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Abstract 



Mh' Unique transformation properties under the hyperspherical inversion of a partial differential 

equation describing a stationary scalar wave in an A''-diniensional (A'^ Jj 2) Maxwell fish-eye 
medium are exploited to construct a closed form of the Green's function for that equation. 

2 ' For those wave numbers for which the Green's function fails to exist, the generalized Green's 

function is derived. Prospective physical applications are mentioned. 
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C^ ■ 1 Introduction 

o ; 

In 1854, Maxwell [1] pointed at a remarkable property of an infinite optical medium with the 
refraction index 

-fe(r) = 4T^ (p>0). (1.1) 

H \ Within the framework of the geometrical optics, he proved that paths of all rays emitted from 

an arbitrarily located point r' are circles having two points in common: the source point r' and 
the image point —p^r'/r'^. For this medium Maxwell coined the name 'the fish-eye'. In 1926, 
Caratheodory [2] observed that there is a geometrical correspondence between the circular rays in 
the fish-eye medium and the geodesies on a sphere. This geometric thread was pursued further by 
other researchers within the group-theoretical framework (cf, e.g., Ref. [3]). 

It seems that Demkov and Ostrovsky [4] were the first to discuss the wavized scalar fish-eye 
problem. Specifically, they considered the equation 



V 



2 , 4zy(iy + l)p2l 



(r2 + p2) 



2\2 



*(r)=0 (p>0) (1.2) 



in R^, subject to the boundary condition that ^(r) vanishes at infinity. They proceeded in two 
directions. First, they solved analytically a spectral problem with v being an eigenparameter 
and showed that the resulting spectrum is purely discrete and eigenfunctions may be expressed 
in terms of the Gcgenbauer polynomials. Second, they proved that Eq. (1.2) possesses a certain 
remarkable transformation property under the geometrical inversion in a certain class of spheres. 



and ingeniously exploited this fact to construct a closed form of the relevant Green's function in M^. 
Group-theoretical properties of the scalar fish-eye wave equation (1.2) were then investigated in R^ 
by Frank et al. [5, 6]. Lately, the two-dimensional fish-eye medium has been studied by Makowski 
and Gorska [7] in the context of the construction of pertinent coherent states. Finally, in two very 
recent papers, Leonhardt [8] and Leonhardt and Philbin [9] have argued that the geometric-optical 
perfect focusing property of the fish-eye medium, discovered by Maxwell, holds as well within the 
wave-optics framework; that issue will be critically reexamined in our upcoming work, with the 
aid of the results presented below. 

The present paper is the first out of a series of several reports in which we shall expose results 
of our research on the wave properties of the Maxwell fish-eye and related media. Here, we derive 
a closed form of the Green's function for the scalar fish-eye wave equation in M.^ , N ^ 2. The 
particular method we employ generalizes the aforementioned one used by Demkov and Ostrovsky 
in the case of TV = 3 and exploits a peculiar transformation property of the TV-dimensional fish-eye 
equation under the hyperspherical inversion. 

The structure of the paper is as follows. In Section 2, we investigate transformation properties 
of a class of partial differential equations under the hyperspherical inversion, with a special focus 
on the TV-dimensional fish-eye equation. The results of that investigation are used in Section 3 
to construct the Green's function for the fish-eye problem. The case when the Green's function 
fails to exist and is to be replaced by the generalized Green's function is considered in Section 
4. Prospective physical applications of the results are briefly discussed in Section 5. The pa- 
per ends with an appendix, in which a number of closed-form representations of the derivative 

[dp-^/^+^{x)/diy]„^n+N/2-i, with X G (-1, 1), A^ G N \ {0, 1} and n G N, required in Section 4, 
are displayed. The list of references attached has been intended to contain representative items 
rather than to be a comprehensive one. An exhaustive listing of works relevant to the Maxwell 
fish-eye problem will be included in one of our forthcoming papers. 

2 Transformation properties of a class of partial differential 
equations under the hyperspherical inversion 

At first, we establish the following 

Lemma 1. If'${r) {with r G M.^ , TV ^ 2) satisfies the equation 

[V^ + fc2n2(r)] *(r) = 0, (2.1) 

then for arbitrary R E K. and a,b E M.^ it holds that 



^ -^(r^^{r-a) + b]=0. (2.2) 



^ + I M '^ I T^{r -a) + b 

Proof. The transformation 

results in the following alteration of the infinitesimal line element: 

Hence, substitution (2.3) implies the following transformation of the TV-dimensional Laplace oper- 
ator: 

v^^t^v.(^^!^v)^ 

Using the easily provable differential identity 

^2{N-2) \ p>2(N~2) I Jl2{N-2) 



' 1^ _ Q^\2{N-2) I l« ^\N-2 ^ l„ „IW-2 l„ „l/V-2 1 ^ l„ „IA^-2 ) ' V^'*^/ 



and exploiting the fact that the function l/|i — a|^ ^is harmonic in K.^, we see that transformation 
(2.3) changes Eq. (2.1) into 



|A'+2 



i?4 



\r — a 



N-2 



J 2 2 

k n 



i?2 



r -a 



:{r-a) + b 



f 



i?2 



r -a 



:{r-a)+b]=0, (2.7) 



which immediately leads to Eq. (2.2). D 

Actually, the above lemma offers a bit more than necessary for the purposes of this paper. In 
view of our needs, in what follows we shall restrict ourselves to the special case when the vectors 
a and b are equal. It is then evident that the resulting transformation 



r n- 



R^ 



;{r - a) + a 



(2.8) 



is the geometric inversion in the hypersphere of radius R centered at the point with the radius 
vector r = a. (In fact, if b = a and fc = 0, the lemma is simply an A^-dimensional extension of the 
well-known Kelvin inversion theorem for harmonic functions [10].) 

Now we turn to the fish-eye problem. Application of the following special case of inversion 

(2.8): 



p2 + a% 

r 1-^ —(r — a) + a => r i-^ 

\r — aM 



r — a\ 



r + a- 



P 



(2.9) 



to the fish-eye refraction index (1.1) gives 



"-fc 



p'- + a^ \ |r — ap 2nop^ 



r -a 



p^ + a^ r^ + p^ p^ + a? 



n[c{r). 



(2.10) 



Combining this property of the index nfo(r) with the result stated in the lemma, we arrive at 
Corollary. // the function '^{r) (with r G R^; N ^ 2) solves the fish-eye equation 

,2 , ^nl^p" " 



(r2 + p'^Y 



*(r-) = {p> 0), 



then for arbitrary a G R the function 



X(a, Vp^+a^)*(r-) 



also solves this equation, i.e., it holds that 



\N-2 



* ( ^ '^{r-a) + a 



(2.11) 



(2.12) 



Anlk'^p^ 
(r^ + p^)^ 



1 



|Ar-2 



^[ .^ ^" (r-a) + a) ^0 (p > 0). (2.13) 



In the particular case of A^ = 3, the above result was established by Demkov and Ostrovsky [4] 
(see also [11]). 



3 The Green's function for the fish-eye problem 

We are now ready to construct the A^-dimensional fish-eye Green's function G,^{r,r'). According 
to the general theory of Green's functions for elliptic partial differential operators, G„{r,r') is a 
single- valued solution to the fish-eye equation 



4 2l2 4 
j2 I ^"0^ P 

(r2 +p2)2 



G,(r,r') = (p>0) 



(3.1) 



everywhere in R except for the source point r = r' , where it diverges according to 



G,(r,r')''-^ —In |r-r'| (iV = 2) 
Zn 



or 



GAr,r') 



1 



In the last equation 



(iV-2)5'Ar_i|r-r'|^-2 
2^^/2 



(A^^3). 



Sn-1 



T{§) 



is a surface area of a unit {N — l)-diniensional sphere S^ ^ embedded in M^. 
require 



GAr,r' 



„N-2 



(3.2) 
(3.3) 

(3.4) 
At infinity, we 

(3.5) 



(the non-zero constant Ci,{r') appearing in condition (3.5) will be determined later). The parameter 
V is defined as 



-l + v/l + 4ngfc2p2 
V = [v ^f \j tor no -^ 0) 



(3.6) 



and reasons for its introduction will become clear shortly. In what follows, we admit that the 
product rL^k"^ 1 hence also v, may be complex. 

At first, consider the case when the source is located at the center of symmetry of the medium. 
Evidently, the corresponding Green's function G'^(r,0) must be spherically symmetric, being a 
function of r = \r\ only. Hence, it follows that Giy{r, 0) obeys 






N-1 d 4i^(j^ + l)p2 



r dr (r^ + p'^y 



G,(r, 0)^0 



(3.7) 



except for the point r = 0, where it behaves according to Eqs. (3.2) or (3.3) with r' = 0. The 

substitution 

.p-. N/2-1 fr'^-p^ 
G.(r,0) = f^ F -^ 



-) 



r^ + p 



(3.8) 



leads to the following differential equation for the function F: 



9n d^ d 

da;^ Ax 1 — x^ 



M 



Fix) = 0, 



where 



and 



X = 



2 2 

r — p 

J.2 _|_ p2 



Ai 



(-1 < a; s^ 1) 



iV-2 



(3.9) 



(3.10) 



(3.11) 



Equation (3.9) is the associated Legendre equation. Its general solution, written in the form 
most suitable for the present purposes, is 



F{x) = AP-^'{x) + BRi^ix) 
with A, B being arbitrary constants and with 



Kix) ^ Q^Jx) + '^Pi^ix) = — ^ 
2 2sm(7r/i) 



e'-P.^(.)-^|^^±^P-(.) 



(3.12) 



(3.13) 



Here, Py{x) and Q'^{x) are the associated Legendre functions (on the cut — 1 ^ a; ^ 1) of the first 
and second kinds, respectively (occasionally, R'^{x) is called the associated Legendre function of 



the third kind). The general character of solution (3.12) follows from the fact that the Wronskian 
of P-''(x) andi?(;(x) is 

WlP-^ixlKix)]^"-^^^, (3.14) 

i.e., it vanishes nowhere. Henceforth, we shall adopt the standard convention and shall write Pi,{x) 
and Ru{x) in place of P^{x) and R^{x). 

With the general solution to Eq. (3.9) in hand, we see that the Green's function Giy{r, 0) is of 
the form 

O.ir. 0) . ,4 (£) ""- p-»/- (^) + B (^) "'" nr-' ('4-4)- (3.15) 

\r/ \r + p / \r / \r ^ p ) 

V being defined in Eq. (3.6). We shall fix values of the constants A and B in two steps. At first, we 
investigate the asymptotics of the expression in Eq. (3.15) as r ^ oo. Using the known formulas 
[12, p. 196] 

P-(.) ^^-" j^ (^) ''' (. ^ -1, -2, . . .), (3.16) 

i?.(x)"^"-iln(l-.T), (3.17) 

K{^Y^''\^'^''n\A(^-^f) ' (ReM>0), (3.18) 

we see that the constraint (3.5) is fulfilled iff i? = 0, and consequently 

G.(.,0) = A (£)""-' P-«>(ll^). ,3.19, 

In the second step, we investigate the asymptotics of the right-hand side of Eq. (3.19) as r ^> 0. 
Exploiting the formulas [12, p. 197] 

P^(,)-_z4+"5i!lMin(l+,) (3.20) 

TT 

and 

^"^ ^"^ ^ WTJ^TWF^7T7}[~^) (^-^>^)' (3.21) 

we find that the constraints (3.2) and (3.3) will be satisfied iff 



1 p^ ( r'-p' 
4sin(7ri^) \r^ + p- 



G.(r-,0) = ^— 73^P. (^-^) (iV = 2) (3.22) 



r(f + .)r(f-.-i)-'' yr^ + p2 



and 

p-N/2+l ( ^ -P 

G.(r,0)^~ ^-^ ";;;/2 " ^' ^ (^rp)^i2-r ' (^^3), (3.23) 

respectively. Since it holds that 

-^"(-)-- r(. + r)r(-.) ' (^-^4) 

Eqs. (3.22) and (3.23) may be collected into a single formula 

p-N/2+l (''__P_ 

G.(r,0) = - ^'^"':l' '^ . .^JJ'' - (3.25) 



4j,N/2 ^^p-jN/ 



From this, using relation (3.16), we deduce that for r' = the constant Ci, in the asymptotic 
relation (3.5) is 

To find the Green's function for an arbitrary location of the source point r', we consider the 
transformed function 

^2 1 „2 



1 ^ /p' + a 



I{a, vV+^)G.(r, 0) = n— G, /^^^(r - a) + a, , (3.27) 



with the center of the inversion sphere (of radius \/ p^ + a^) located at the point 

a = -r'l-^. (3.28) 

Using Eq. (3.25), the explicit form of this transformed function, denoted hereafter as g^(r,r'), is 
seen to be 

p-N/2+i{ , , 2p\r-r'f \ 



47r^/2 Vp^ 



lr-r'r/2-i 



^ + ^72 



Af/2-l 



(3.29) 

In view of the results of Section 2, we know for sure that the function in Eq. (3.29) solves the 
fish-eye equation, except, possibly, for some isolated points. It is evident that the points at which 
the behavior oi g^{r,r') should be investigated are the two finite points r — r' and r — ~r' p^ /r'^ , 
and also the point at infinity. Using the asymptotic relations (3.20) and (3.21), we derive 

5,(r,r')''-^';^ln|r-r'| (TV = 2) (3.30) 

and 

1 / r' \ ^"2 

"-'-■-'' --^- (iv-2)5.-.i.-.r-4 ^) '"*''• <"" 

i.e., the 'inverted' function diverges for r — > r' in the same manner (save for the factor (r'/p^)^^^ 
when iV ^ 3) as the Green's function Gu{r,r') (cf. Eqs. (3.2) and (3.3)). Furthermore, it is seen 
that for r — > 00 the function (3.29) decays asymptotically as 

p-Ar/2+l f P ^ 



gAry)-—- -- — /, ^(^^J ,.^2 ' ^ (3-32) 

i.e., in the same functional manner with r as prescribed for Giy{r,r') in Eq. (3.5). Finally, with 
the help of the identity 

/ 2 \ 2 

„ 9/ ,x9 2r'2 ir + r'-^] 



(r2 +p2)(j,'2 _|_p2-) (^2 ^p2-)(^/2 _|_p2-) 

and the asymptotic relation (3.16), we find that for r -^ —r'p^/r'"^ the function (3.29) remains 
finite: 



r.r^-r'A/r- F (f + .) F (f - . - l) // ^ t"'^) 1 

4^A^/2r(f) \p) (^'2+^2)Ar-2- 



9Ary) -^' '' ._:n.y^N. -{-) TZPT-^KW^,- (3-34) 



Thus, we see that the function {p'^ /r')^^'^g^{r,r') satisfies all conditions imposed on G^{r,r') in 
Eqs. (3.1)-(3.5). Hence, we conclude that the closed form of the TV-dimensional fish-eye Green's 
function for an arbitrary location of the source point r' is 

p-Nl2+i( . I 2p^{r-r'f 
0^(,r,rj- —TTjT — [6.6i)) 



4t,n/2 



\r ~ r 



'|W/2-i ^ I ^'i^ 



and that the constant C^{r') in the asymptotic constraint (3.5) is 

c(r')- r(f + ^)r(f-.-i) .px^/2-i 




N/2-1 



(3.36) 



Since the differential operator in Eq. (3.1) is symmetric with respect to the scalar product 
(x 0)/v = /r" d^r x(r)(/)(r), the fish-eye Green's function should be symmetric with respect to 
the interchange of the source and observation points: 



G,(r,r') = G.(r',r). 



(3.37) 



However, the representation of G„{r, r') given in Eq. (3.35) does not exhibit this property explicitly. 
To show that nevertheless relation (3.37) is satisfied, we observe that it holds that 



,P 



Vr^r'^ + 2p^r ■ r' + p'^. 



(3.38) 



Consequently, the Green's function (3.35) may be alternatively rewritten in either of the following 
two manifestly symmetric forms: 



G,{ry) 



N/2-1 p-N/2+l f -I 



2p^(r -r 



l\2 



(r2 +p2)(y,'2 ^p2-) 



47r^/2 



|r-r'|"/2-i(rV2 + 2p2r-r' + p4) 



W/4-1/2 



(3.39) 



or 



G.(r,r') 



r(f + .)r(f-.-i 

47r^/2 



pN/2-lp-N/2+l(^^ 



2p\r-r'f 

(r2 +p2)(^/2 _j_ p2-^ 



(r/)JV/4-l/2|^_^/|W/2-l 



N/4-1/2 



.P 



N/4-1/2' 



(3.40) 



Still another representation of Gi,{r, r') displaying its symmetry is the one in terms of the Gegen- 
bauer function of the first kind. Using the known relationship 

^"(^) - 2&j^m^)^' " .^)-/2.V4p-^..V2^(,) (_i ^ , ^ 1), (3.41) 

Eq. (3.35) is transformed into 



■)N-4r ( N-l 



G^,{r,r') 



r(^) 



P ^v-N/2+1 I ^ 



2p'{r-r'f 

(r2 +p2)(^/2 ^p2-) 



T,(N-l)/2 sin [tt (^ - v)] (r2 + p2)Ar/2-l(^,2 + ^2)Ar/2-l 



(3.42) 



Let us consider some particular cases. For N = 2, from cither of Eqs. (3.35), (3.39) or (3.40), 
with the aid of Eq. (3.24), we find 



GAr,r') 



1 



4sin(7ri^) 



P. 



2p'{r~r'f 

(j,2 _^ p2^(j.l2 _|_p2) 



{N = 2). 



(3.43) 



Next, it appears that the representations of G^{r, r') found above simphfy greatly when TV is odd, 
as then the Legendre function appearing in Eqs. (3.35), (3.39) and (3.40) may be expressed in 
terms of trigonometric and inverse trigonometric functions either as [12, p. 168] 



p-'^/'^Hx) 



or as [12, p. 169] 



l^-j' (. 2N-W/4+1/2 V- . xfc r(fc + t/-- + 2) 

2^/2-2^1 ■ ) 2^^ y > kW{k + iy+l){^-k)l 



fc=0 



sin [(2fc + i^ - f + 2) arccos x] {N odd, iV > 3) 



(3.44) 



(Af-3)/2 



p-^/^+i(.) = ^^r(.-f+2) y: 

sin [(fc 



(fc + ^)! 



..0 2mT{k + .+ l){^^ky. 

{N odd, N ^3). 



i^ + i) arccos X + (fc 



^)f] 



(^1 _ 2,2-)fc/2+l/4 



In the simplest case of A^ = 3, we have 



P-'/'{:c) = J±{l-x') 



o.-i/4sin[(j/+i) 



arccos x 



(3.45) 



(3.46) 



Using this representation of Pi, {x) in Eq. (3.39), the latter being specialized to the case iV = 3, 
after some straightforward movements and with the help of the identity 



r(^ + ^)r(i-z^) 



cos(7ri^) ' 



(3.47) 



we arrive at 



GAr,r') = - 



Vp + /?)P^ + /?) 



47rcos(7ri^) \r - r'\^Jr'^r''^ + Ip^r ■ r' + p'^ 



X sm 



1\ / 2p\r-r'f 
ly + - ] arccos -1 + —^ „,, ,„ ^ 



(TV = 3). (3.48) 



Equivalence between the result in Eq. (3.48) and the following expression (modified to conform 
with the present notation and corrected for a sign error): 



GAr,r') = 



1 



VP + p^)P^ + /^ 



47rcos(7rz^) \r — r'\^r'^r''^ + 2p^r ■ r' + p"^ 
y'r^r'^ + 2p2r • r' + p^ 



X sm 



(2i/ +1) arctan- 



p\r — r'\ 



{N ^ 3), (3.49) 



given by Demkov and Ostrovsky in Ref. [4] , may be easily established with the aid of the well-known 
inverse trigonometric identity 



1 1-^2 
arctan^ — - arccos — (C 5^ 0)- 



(3.50) 



4 The generalized Green's function for the fish-eye problem 

A glance at either of Eqs. (3.35), (3.39) or (3.40) reveals that the fish-eye Green's function G,^{r, r') 
fails to exist for the following values of i/: 



V ^n + N/2 - 1 



or V — —n 



N/2 (n e N), 



(4.1) 



which, by the way, are solutions of the quadratic equation 



1^(1^ + 1) = [n + YJ (" + y-i 



(4.2) 



If either of the conditions set in Eq. (4.1) holds, one seeks the generalized Green's function 
G'„+jv/2-i ('":'"') = G_n_j^/2{r,r'), defined through the limiting relation 



Gn+N/2-i{r,r') 



lini 



i/(!/+l)-!-(n+Af/2)(Ti+Af/2-l) 



d 



d[iy{iy + l)] 
1 



K^ + 1)- {n+-j ("+Y-1 



lim 



d 



2n + N — 1 u^n+N/2-i dv 



G,{r,r') 
(iy-n- — + lj L + n+—jG^{r,r') 



(4.3) 



If, for instance, representation (3.39) of G^{r,r') is used in Eq. (4.3), this results in 



G„+^,/2_i(r,r') = (-)' 



(n + iV-2)! 



JV/2-1 



47r^/2n! 



ip{n + N ~1) - tp{n + 1) + 



1 



|r-r'|^/2-i(r2r'2 + 2p2r.r' + p4) 



2p'(r-r'f 



N/i-l/2 



Qp-N/2+ll_^ 



2n + N -I 

2l„ „l\2 



p-N/2+1 , _, 
^n+N/2-1 ' ^ 



(r2 +p2)(^/2 ^^2) 



2p^{r~r'y 

(j,2 _|_^2-)(-^/2 ^_p2-) 



dv 



where 



V'(C) 



u=n+N/2-l, 

1 dr(c) 



(4.4) 



(4.5) 



r(C) dc 

is the digamma function. A number of closed-form representations of the derivative 
[dPi, {x) / di>\i,^n+fq /2-1 required in Eq. (4.4) may be derived from the author's findings for 

[dP^"^{z)/di']^=n, ^ e C\ (—1,1), presented in Refs. [13, 14]; the simplest, and thus potentially 
most useful, of these expressions are listed in Appendix A. 
In the particular case of A^ = 3, Eq. (4.4) yields simply 



G„+i/2(r,r') 



(-)" V(^'+P')(^"+P') 



47r2 |r -r'|^r2r'2 +2p'^r ■ r' + p^ 

( 2p^{r-r'f 

{n + 1) arccos ( -1 + ,,2^ ow ,2^.2 



sm 



(n + 1) arccos — 1 



(r2 +p2)(^'2 _|_p2-) 

2p'{r-r'f 

(r2 +p2)(j./2 _|_p2) 



arccos — 1 



2p%r-r') 



l\2 



(r2 +p2)(^/2 +p2-)^ 

(AT = 3). (4.6) 



2(n+l) J 

5 Prospective applications 

The closed-form representations of the fish-eye Green's function G,^{r,r') and of its generalized 
counterpart G„^N/2-i{r,r'), found in this work, are certainly interesting for their own mathe- 
matical sake. They appear, however, to be also useful in the physical context. In a forthcoming 
report, we shall use them to show that, despite of claims to the contrary [4, 8, 9], in wave optics 
the infinite Maxwell fish-eye medium does not possess the same perfect focusing properties as it 
has in geometrical optics. Next, it has been confirmed [15] that the use of either of the closed-form 



expressions for G^{r,r') listed in Section 3 simplifies greatly the mathematical analysis of wave- 
optical properties of cylindrical (A^ — 2) and spherical (A^ — 3) gradient-index lenses with the 
fish-eye refraction index (1.1) and of finite radii Hens ^ p\/2rio ~ 1- Finally, in yet another forth- 
coming paper, we shall show that there is a close mathematical relationship between the wavized 
Maxwell fish-eye problem in M^ and the A'^-dimensional Schrodinger-Coulomb problem in momen- 
tum space; in particular, we shall provide there an integral expression for the momentum-space 
Schrodinger-Coulomb Green's function in terms of the fish-eye Green's function discussed above. 



A Appendix: The derivatives [dP,, 
iVGN\{0,l} 



-N/2+1 



x)/diy]j,=n+N/2-i for 



From the relations 



and 



P-"(X) ^ c-'™/2i!i-^P™(x + iO) (0 ^ m ^ n) 



(A.1) 



diy 



-i7rm/2 



(n-m)!9P^(a;-hiO) 



(n + my. dv 



[il^{n + m+l)-ip{n-m+ l)]p-"(x) 
(0 ^ m :^ n), (A.2) 



and from a number of closed-form expressions for the derivative [dP^ {z) / dv\^=n, with 2 S C \ 
(—1,1) and ^ 771 ^ n, found by the present author in Refs. [13, 14], one may derive, among 

others, the following representations of [dPu ' {x) / dv\ij^n+N /2-1, with x e [—1, 1] and with N 
being an even natural number greater than zero'} 



dP, 



-N/2+1 



(x) 



dv 



^ P 



-N/2+1 



l-\-x 



-N/2+1 



iy=n+N/2-l 



+ 



n\ 



{n + N -2) 



\:N/2-li^) In ^ - 2V^(^ + ^ - l)^nT^2-l(-) 



l-x^^"/^-^/^^^_^, (fc + n + iV-2)! 



fc=0 



fc!(fc + iV/2-l)!(n-fc)! 



X [27p{k + n + N -l)-ii}{k + N/2)] 



1-x 



^_^^N/4-i/2n+m-i ik-^-n + N/2-imk + N/2) 

^ ^ ' k\{k-^N/2-iy.{n + N/2-k-iy. 



1 + x 



fe=0 



l~x 



(A.3) 



dP, 



-N/2+1 



(x) 



P 



-N/2+1 



Af/2-lV^)^° 



1+.T 



Af/2+1 



i^=n+N/2-l 

n 



I /l-x2^^/^~^/^ " 



{n + N -2y 



E(-) 



fe=0 



1-x 

1 + x 



N/4-1/2 n+N/2-1 



E (-) 



2i,{n + N/2)p-l'^^;^^{x) 

f^{k + n + N -2y.tP{k + N/2) (1-x 
k\{k + N/2-iy.{n-ky. 

{k + n + N/2-iy. 



k=0 



k\{k + N/2 - l)!(n + N/2 - fc - 1)! 



X [2ilj{k + n + N/2) - iP{k + N/2)] 



1-x 



(A.4) 



^Attention! The Mathematica 7.0.0 function Legendr eP [nu , mu , 2 , x] evaluates incorrectly (the modulus is correct 
but the sign is wrong) numerical values of the associated Lcgcndre functions ^2~ri+l~ (^) ^^'^ ^2n+2 (^) ^°^ « € N, 
— 1 ^ a; ^ 1. An empirically discovered remedy is to subject the variable x to the action of the function SetPrecision 
before it is used as an argument of LegendreP. 
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ap-^/'+^(x) 



dv 



_ p-N/2+l ,. , 1^ 



v=n+N/2-l 



-N/2+1 



[i,{n + iV - 1) + ^(n + 7V/2)]P„7^;;_\(x) 



(n + iV-2)! 



^_^2\ N/i-l/2 n 



E(- 



fc=0 



fe (A: + n + TV - 2)!V'(A: + n + 7V-l) /^l-x 
' A:!(fc + iV/2-l)!(n-A:)! 



l-^^^^''^^'"^v^'''( ^, (fc + n+^JV/2-l)!^(fc + n + JV/2)^ f ^'^^ V, (A.5) 



1+.T 



fe=0 



A:!(fc + N/2 - l)!(n + Ar/2 - A: - 1)! 



dP, 



-N/2+l 



(x) 



9:^ 



= P 



-N/2+l 



n+N/2- 



^{x) In ■ 



l + x 



i/=n+N/2-l 



( ^„jl-x^y"^^"^/^"g^fc + n)!(7V/2-fc-2)!^l + . 



*;=0 



fc!(n + Ar-fc-2)! 



+ (-)" 



1 + a; 
l-x 



Af/4-1/2 "+W/2-1 



E (- 



(A: + n + 7V/2-l)! 



*;=0 



fc!(fc + 7V/2 - l)!(n + iV/2 - A: - 1)! 
fe 



X [2?A(fc + n + Ar/2) - V(fc + iV/2) - V'lfc + 1)] 



1 + x 



(A.6) 



ap-"/^+\x 



dv 



p-^/^+' (x) In i±^ 



,y=ri+Ar/2-l 



[^(„ + TV - 1) - ^(n + l)]p-+ji't\{x) 

n! /l-a;\'^/*"'/'T^'(fc + n + 7V/2-l)!(iV/2-/fc-2)! /l + a;\'= 



-(-) 



(-)" 



(n + iV-2)! V1 + .T 
' /1-a; 



{n + N -2)1 



E 

2 \ Af/4-1/2 n 



E(- 

fc=0 



A:!(n + iV/2- fc-l)! 

fc (fc + n + JV-2)! 
A:!(fc + Ar/2-l)!(n-/c)! 



X [2i/'(fc + n + iV - 1) - ^(/c + 7V/2) - -(/-(fc + 1)] 



1 + x 



(A.7) 



9p-^/^+^(x 



9z/ 



u=n+N/2-l 



p-N/2+1 U)l^l±^ 



,-N/2+l 



l + a;\ / 1 — x^ 



l-x, 
(fc-1)! 



l-x 



+ [^(n + 1) + VXn + A^/2)]P„-_;^;;;_\(x) 
-{-Y'n\{n + N/2-l)\ 

N/2-l 

^ (fc + n)!(fc + n + Ar/2 - l)\{N/2 - k - 1)\ \l + x ^ 
A fe V(^ + A^/2- A:) +1/^(71- fc+1) fl-x^^ 

fc=0 



k\(k + 7V/2 - l)!(n - k)\{n + 7V/2 - /c - 1)! \l + x 



(A., 
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ap-^/'+^(x) 



p 



-Ar/2+l 



u=n+N/2-l 



W/2-l(^)l'^' 



,-Af/2+l 



1 + a: 



+ [V^(n + 1) + VXn + 7v/2)]p„7;;;_\(x) 

-(-)"n!(n + A^/2-l)! 

Af/2-2 

X E 

fe=0 



T X Af/4-1/2 /, \ n+N/2-1 



l + x^ 
{N/2-k-2)\ 



l + x 
k\{n + N/2 -k- l)!(n + A^ - fc - 2)! \l - x 



{-Yn\{n + N/2-l)\ 

n 



, , X Af/4-1/2 /, \ n+W/2-1 



^(fc + l)+^(fc + Ar/2) 



fc=0 



l + x 
kl{k + N/2 - l)!(n - ky.{n + N/2 - k - 1)1 ^1 - x 



(A.9) 



ap-^/'+\x) 



9j^ 



P 



-Af/2+l 



n+N/2- 



^(x)ln 



l + x 



;/=n+Af/2-l 



+ [2^(2n + N-l)-i,in + N-l)-^in + N/2)]p-^^^+\ix) 



EK 



k-\-7i 



2k + N -1 



fe=0 
Ar/2-2 
y^ ^_\k+n+N/2 

fc=0 



(n-fc)(fc + n + AT- 1) 



n!(fc + Ar-2)! 



/c!(n + Ar-2)! 



P 



-Af/2+l 



k+N/2- 



\{x) 



2k + 1 



-P, 



-N/2+1 



(x). 



(A.IO) 



(n + Ar/2-/c-l)(fc + n + Ar/2) '' 
If, in turn, A^ is an odd natural number greater than 1, then from Eqs. (3.44) and (3.45) one obtains 



9p-^/^+^(x) 



dh' 



v=n+N/2-l 

H-x^Y 



~^n+N/2-i^"'"' ^-rceosx 



(N~3\, , , (W-3)/2 

_OZi_ h „2\-A^/4+l/2 >^ ^ ^fe 
2^/2-2 /^ I 

V fc=0 



(fc + n)! 



^ k\{k + n+^)\{^-k)\ 
X [V' (fc + n + ■^^) - V'(fc + n + 1)] sin[(2fc + n + 1) arccosx] (A. 11) 



and 



5P-^/'+\x) 



-Qn+iV/t-i(2^) arecosx + V'(ri + l)^„^^/:]l^(2;) 



,y=ri+Ar/2-l 

2- <"f ^/^ (fc+i^)!^(fc + n+^) 
TT ^ 2'=fc!(fc + n+ V)!(i^-fc)! 

sin [(fc + n + ^) arccosx + (fc - ^) f ] 



n - 2,2\fc/2+l/4 



(A.12) 



where 



p-^'^+^ (x) - 

^n+N/2-l\-^l 



\^2-)- (-, _ 2N-^^/4+l/2 V- 
2W/2-2^ {' "^ ) 2^ 

X sin[(2fc + n + 1) arccosx] 



iN-3)/2 



(fc + n)! 



(A.13) 



12 



or, equivalently, 



P 



-N/2+1 



n+N/2- 



l(^) 



(Ar-3)/2 

E 



i\ 



(k + ^y. 



^^ 2^kl{k + n+^)l{^-k)l 

sin [{k + n+ ^) arccosx + (fc - ^) f ] 
n _ 2;2)fe/2+i/4 



(A.14) 



and also 






^-N/2+l , . 
^n+N/2-l\^> 


2^/2-1 ^' ^ ^ ^ ^ ^ fc!(fc + n+\i)!(V- 


-fc)! 




X cos[(2fc + n + 1) arccosa;] 


(A.15) 


or, equivalently, 






r- (JV-3)/2 .,^ ^_3., 






cos [(fc + n + ^) arccosa; + (fc - ^) f " 

^ (1 - a;2)fe/2+l/4 


(A.16) 
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